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ABSTRACT 

Aims. We present an empirical model of Type la supernovae spectro-photometric evolution with time. 

Methods. The model is built using a large data set including light-curves and spectra of both nearby and distant supernovae, the latter 
being observed by the SNLS collaboration. We derive the average spectral sequence of Type la supernovae and their main variability 
components including a color variation law. The model allows us to measure distance moduli in the spectral range 2500 - 8000 A 
with calculable uncertainties, including those arising from variability of spectral features. 

Results. Thanks to the use of high-redshift SNe to model the rest-frame UV spectral energy distribution, we are able to derive 
improved distance estimates for SNe la in the redshift range 0.8 < z < 1.1. The model can also be used to improve spectroscopic 
identification algorithms, and derive photometric redshifts of distant Type la supernovae. 

Key words, supernovae: general - cosmology: observations 



1. Introduction 

The evolution of luminosity or angular distance with redshift is 
an essential observable to constrain the equation of state of dark 
energy, responsible for the acceleration of the expansion of the 
Universe. Type la supernovae (SNe la) are, today, among the 
best distance indicators. They can be observed up to high red- 
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* Based on observations obtained with MegaPrime/MegaCam, a joint 
project of CFHT and CEA/DAPNIA, at the Canada-France-Hawaii 
Telescope (CFHT) which is operated by the National Research Council 
(NRC) of Canada, the Institut National des Sciences de l'Univers of 
the Centre National de la Recherche Scientifique (CNRS) of France, 
and the University of Hawaii. This work is based in part on data 
products produced at the Canadian Astronomy Data Centre as part 
of the Canada-France-Hawaii Telescope Legacy Survey, a collabora- 
tive project of NRC and CNRS. Based on observations obtained at the 
European Southern Observatory using the Very Large Telescope on the 
Cerro Paranal (ESO Large Programme 171.A-0486). Based on observa- 
tions (programs GN-2004A-Q-19, GS-2004A-Q-11, GN-2003B-Q-9, 
and GS-2003B-Q-8) obtained at the Gemini Observatory, which is oper- 
ated by the Association of Universities for Research in Astronomy, Inc., 
under a cooperative agreement with the NSF on behalf of the Gemini 
partnership: the National Science Foundation (United States), the 
Particle Physics and Astronomy Research Council (United Kingdom), 
the National Research Council (Canada), CONICYT (Chile), the 
Australian Research Council (Australia), CNPq (Brazil) and CONICET 
(Argentina). 



shifts, and provide a precise estimator, with a typical dispersion 
of 7% on distan ce, when not limite d by the measurement uncer- 
tainties (see e.g. lAstier et al.ll2006l hereafter A06). 

As the number of SNe la in Hubble diagrams increases, 
systematic uncertainties are becoming the main limitation to 
the accuracy of measurements of cosmological parameters with 
SNe la. Among the potentially serious systematic uncertainties, 
the dominant ones are a possible evolution of the supernova pop- 
ulation, the photometric calibration, the modeling of the instru- 
ment response and the uncertainties arising from SN la large 
spectral features, including their possible supernova to super- 
nova variations. In this paper, we aim at addressing principally 
the latter although the proposed method makes it also easy to 
account for the modeling of the instrument and to propagate the 
model uncertainties. 

Various approaches to distance estimation have been pro- 
posed, using l ight-curve shape parameters (Amis or a stretch 



posed, using light-curve shape parameters (Amis or a str etcf 
factor , see e.g. lPhillipsl [i993; Rie ss et al.l[l9 95; Perlmutte r et al 



I 1997I) or colo r information (TWang et all l2003bi IWang et al . 
I2005I) . or both dRiess et al.lll996l lTripplll998l iGuv et al.ll2005l) . 
None of these methods really address the problem of uncertain- 
ties due to the variability of the large features of SNe la spec- 
tra. In this r espect, most method s rely on the spectral sequence 
provided bv lNugent et al.l d2002l) (hereafter N02). This is a seri- 
ous concern since this could possibly result in sizable (common) 
system atic effects in the distance measurements. In IGuv et all 
(2005) (hereafter SALT), we have applied broadband corrections 
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to the spectral sequence N02 as a function of phase, wavelength 
and a stretch factor so that the spectra integrated in response 
functions match the observed light-curves. This provided us with 
a tool to fit the observed light-curves without correcting the data 
points since the /if-corrections were naturally built into the model 
(see e.g. N02 for a definition of /f-corrections). This approach 
has been quite successful when applied to estimating SNe la dis- 
tances at high redshift. However it did not address the problem 
of variability of spectral features either. 

In this paper, we use a similar framework. Whereas we only 
used multi-band light curves to train the model in SALT, here we 
include spectroscopic data to improve the model resolution in 
wavelength space, be fully independent of the spectral sequence 
N02, and more generally extract the maximum amount of in- 
formation from the current data sets. Modeling the supernova 
signal in spectroscopic space ensures that the /f-corrections are 
treated in a consistent manner since there is a single model to 
address both light-curves and spectra. It also permits a coherent 
propagation of errors, from the fit of light-curves to distance es- 
timate. The model is allowed to vary as a function of phase and 
wavelength with a small number of a priori unknown intrinsic 
parameters and a color variation law which is also adjusted dur- 
ing the training process. The main goal of this approach is to 
provide the best "average" spectral sequence and the principal 
components responsible for the diversity of SNe la, so that the 
model can account for possible variations in SNe la spectra at 
any given phase. 

The flux normalization of each supernova is a free parameter 
of the model. Hence, we do not need to know their distances to 
train the model. This allows us to use both nearby SNe which are 
not in the Hubble flow and high-redshift ones without any prior 
on cosmology. Using high-redshift supernovae permits to model 
the rest-frame UV emission which is invaluable to improve dis- 
tance estimates of supernovae found at redshifts larger than 0.8. 

In Sect. [2] we present the model implementation. The su- 
pernova data sets used for training the model are described in 
Sect. [3] Some technical aspects of the training procedure are 
then given in Sect. |4j and in Sect. [5] we present some qualita- 
tive aspects of the resulting model. In an effort to improve dis- 
tance estimates for a cosmology application to SNe la surveys, 
we quantify the remaining variability beyond the principal com- 
ponents extracted in Sect. [6] We show how distance estimates 
of SNLS distant SNe are improved with this approach in Sect. [7] 
and discuss several other possible use of the model in Sect. [8] 

2. The Type la supernova spectral sequence model 

We aim at modeling the mean evolution of the spectral energy 
distribution (SED) sequence of SNe la and its variation with a 
few dominant components, including a time independent varia- 
tion with color, whether it is intrinsic or due to extinction by dust 
in the host galaxy (or both). The following functional form for 
the flux is used 

F(SN,p,A) = xo x [Mo(p,A) + x\M\(p,A) + ...] 

x exp[cCL(^)] (1) 

where p is the rest-frame time since the date of maximum lu- 
minosity in B-band (the phase), and A the wavelength in the 
rest-frame of the SN. Mo(p, A) is the average spectral sequence 
whereas M^(p, A), for k > 0, are additional components that de- 
scribe the main variability of SNe la. CL(A) represents the av- 
erage color correction law. As for SALT, the optical depth is 
expressed using a color offset with respect to the average at the 



date maximum luminosity in B-band, c — (B — V)max — {B— V). 
This parametrization models the part of the color variation that 
is independent of phase, whereas the remaining color variation 
with phase is accounted for by the linear components, xq is the 
normalization of the SED sequence, and x^ for k > 0, are the in- 
trinsic parameters of this SN (such as a stretch factor). To sum- 
marize, whereas (M*) and CL are properties of the global model, 
(xk) and c are parameters of a given supernova and hence differ 
from one SN to another. 

Except for the color exponential term, Eq. Q] is equivalent 
to a principal component decomposition. However, a princi- 
pal component analysis cannot be used since this would re- 
quire having an homogeneous and dense set of observations for 
each SN, namely one spectro-photometric spectrum every 4-5 
days, which is not presently available (n ote that current ongo- 
ing SN programs such as the SN factorv. lAldering et all [20021 
the Carnegie Supernova Program. lHamuv et alj|200q, the CfA 
Supernova programQ and the LOTOSS projectQ, should provide 
such data in the coming years). So we resorted to using a method 
able to deal with missing data. The method used is described in 
the next section. 

2.1. Model implementation 

The phase space that we want to model (wavelength range times 
phase range) is not covered by the set of observations of any 
given supernova. We typically have for each supernova a lim- 
ited set of light-curves points observed with different filters and, 
for some supernovae, one or several spectra at different phases. 
However, when using an ensemble of SNe, this phase space can 
be correctly sampled and if the data set is large enough, several 
components can be extracted. 

In order to link the model defined by a limited set of param- 
eters and the SNe observations, we used a basis of functions, 
as function of phase and wavelength [fi(p,A)]. We used third 
order B-splines (to ensure continuous second derivatives). The 
actual choice of the basis is irrelevant in the phase space regions 
which are densely covered by data, as long as it provides a suffi- 
cient resolution to follow the observed variability of the SED se- 
quence as a function of phase and wavelength. Choosing another 
basis will modify the model in regions where it is poorly con- 
strained, such as very early spectra (p < -15 days). As described 
in section [6] those poorly constrained phase space regions are 
identified after the training using a jack-knife technique. In this 
framework, a model is a linear combination of the basis func- 
tions and can be described by a vector M. Each measurement at 
a given phase and wavelength m(p m , A m ) is then compared to the 
model with a vector H„, (with values H m j = fi(p m , A m )), so that 
the expected value for the model at (p m , A m ) is the scalar product 

2.2. The use of spectral information 

Most spectra of SNe la available in the literature are not cali- 
brated photometrically. Their flux calibration have broad-band 
systematic uncertainties. 

One way to circumvent this difficulty consists in photometri- 
cally "re-calibrating" a given spectrum using the available light- 
curves for this particular SN. However since the full SED se- 

1 CfA Supernova Group: 

cfa-www. har\'ard. edu/oir/Research/supernovafindex. html 

2 The Lick Observatory and Tenagra Observatory Supernova 
Searches: astro.berkeley.edu/ baitflotoss.html 
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quence model is needed to derive an accurate interpolation be- 
tween light curve points at the date of the spectroscopic observa- 
tion, and the spectra are needed as well to accurately model the 
spectral features of SNe, the photometric "re-calibration" of the 
spectra has to be included in the global minimization procedure. 
We have chosen to parameterize the "re-calibration" function 
with the exponential of a polynomial (to force positive corrected 
fluxes), with the degree of the polynomial limited by the number 
of light curves for the SN, and the wavelength range of the spec- 
trum. This re-calibration function is applied to the model, for 
which the SN parameters are mostly constrained by the light- 
curves. Thanks to the simultaneous use of a large amount of 
SNe data, we do not need to have photometric observations at 
the same epoch as spectroscopic ones. 

Statistical errors are rarely provided with the spectra. We 
have evaluated them using the fact that all SNe spectral features 
are broadened due to the kinematics of the ejected matter. So we 
expect an intrinsic correlation length greater than 30A (for a ve- 
locity range larger than about 2000 km. s~ 1 ) which permits one to 
evaluate the photon noise in spectra (assumed to be white noise). 
Nonetheless, we scaled errors so that the weight of spectra was 
of order of that of light curves for which we expect lower sys- 
tematic errors. This weighting, along with the resolution of the 
re-calibration function, is a bit arbitrary but can not be avoided 
at this stage due to the quality of currently available data sets. 

3. The training supernova data sets 

In this section we describe the data sets used for training the 
model. 

In the proposed model (Eq. [TJ, the overall SED sequence 
normalization of a given SN is a free parameter (xo). As a conse- 
quence, it is possible to use both the very nearby supernovae data 
that are not in the Hubble flow (z < 0.001) and the high redshift 
ones without adopting values of cosmological parameters. 

Nearby supernovae have much higher signal to noise than 
their distant counterparts over a much wider range of phases. 
One important difficulty however in using nearby SNe is that 
they suffer from potentially large systematic errors in the ultra- 
violet (UV), since the atmospheric extinction is strong and dif- 
ficult to model in this wavelength range. Including SNe from 
a large redshift range helps to sample homogeneously the rest- 
frame visible wavelength range with both photometric and spec- 
troscopic data, especially in the rest-frame UV. Indeed, if only 
nearby SNe are used, photometric data do not cover the gaps 
between the central wavelength of the filter set used (mostly 
Johnson-Cousins UBVRI), and one relies only on spectra for in- 
terpolation between those bands, which may introduce (weak) 
systematic effects on -corrections. 

One may argue that possible evolution of the SNe la with 
redshift might cause some problems with the modeling since 
objects at all redshifts are used to obtain the model. Actually, 
the model describes an average SN la at an average redshift but 
evolution can still be studied. For instance, without any a priori 
on the effect of evolution on SED sequence, one can look at the 
variation of the (xi) parameters with redshift. 

3. 1 . The nearby supernova sample 

We use a sample of 52 nearby supernovae (without restricting 
ourselves to very nearby ones as in SALT) listed table [2] Those 
SNe were selected from the quality of their light-curve sampling 
where we basically require measurements before the date of 



maximum to ensure a good estimate of the luminos ity at peak. A 
large fr ac tion of thos e SNe l ight- curves come fromlHamuy et al.l 
(Il996bl) : iRiess et alj (Il996l) and iJha et al.l d2006l) . We did not 
consider 1991bg-like SNe la (with very low stretches). They 
have such different light-curves and spectra that the linear model 
we consider cannot fit those along with other SNe la. This is not 
a problem since we aim at modeling the bulk of the SNe la pop- 
ulation (and we do not expect to detect many of those objects at 
high redshift). 

We do not use any spectra without photometric data for the 
same SN (at least two light-curves in different filters), so that the 
date of maximum, color and (xk) can be determined. However, 
since spectra are calibrated on the model and not on the pho- 
tometric data, we do not need simultaneous photometric obser- 
vations; we just need enough photometric observations to de- 
rive the SN parameters. From the sample of 52 SNe, we were 
able to gather 264 spectra for 16 SNe. There are 10 spectral se- 
quences (with more than 10 spectra), namely 1989B, 1990N, 
1991T, 1992A, 1994D, 1996X, 1998aq, 1998bu, 1999ee, and 
2002bo (see Table [2] for the complete list of spectra and their 
references). 

All a vailable UV spectra from the International Ultraviolet 
Explorer dlNESll2006l) were included. This is very helpful since 
most high-redshift SNe spectra which cover the rest-frame UV 
range have a low signal to noise ratio. For all spectra from 
ground-based observations, we do not consider any measure- 
ment below 3400 A because of the strength of the atmospheric 
absorption in this spectral region. 

3.2. The high-redshift supernova sample 

We used a set of 121 Type la supernovae light-curves obtained 
by the Supernova Legacy Survey (SNLS) during the first 2 years 
of the survey (see Table|2]). The light curves were obtained with 
the same reduction pipeline as described in A06, but with new 
images in the photometric fit, so that light-curves have more data 
points and with improved statistical accuracy since the reference 
data used to anchor the estimate of the galaxy brightness be- 
come deeper with time (thanks to the rolling search observing 
strategy). All the 71 SNe used for cosmology in A06 were used, 
with 50 additional ones. 

In addition to the light-curve point s, we used 39 h igh- 
redshift SNLS spectra obta ined at VLT dBasa et al.ll2007l) and 
Gemini dHowell et alj|2005h during the regular SNLS observa- 
tion programs, which aim at typing and measuring the redshifts 
of SN candidates. Obviously more spectra were recorded (at 
least one for each SN) but we choose to use only those with neg- 
ligible residual contamination from the host galaxy. The contri- 
bution from the host galax y was removed i n the reduction proce- 
dure of the VLT spectra dBaumontll2007l) . For all spectra, the 
remaining contamination was evaluated a posteriori using the 
model itself with the following procedure: using the SN parame- 
ters retrieved from the fit of light-curves, the model was fit on the 
SN spectrum with re-calibration parameters, and an additional 
contribution of the host galaxy. We used for this purpose tem- 
plates of elliptical, SO, Sa, Sb and Sc galaxies, the actual galaxy 
type was fitted at the same time as its normalization. All spectra 
with a non zero contribution of the galaxy (at 68% confidence 
level) were not used in the training sample. 

Figure Q] shows the (p, A) phase-space region covered by the 
photometric and spectroscopic data sets. Since we do not use 
infra-red photometric data, the re-calibration of spectra may not 
be reliable for rest-frame wavelengths larger than 8000 A, which 
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Fig. 1 Phase-space mapping by photometric data (top) and spec- 
tra (bottom). For photometric observations, the rest-frame cen- 
tral wavelength of the filter is considered. 



is the central wavelength of the /-band filter. Also, we have lit- 
tle spectroscopic information in the UV for phases earlier than 
-10 days or greater than 10 days since the spectroscopic obser- 
vations of the SNLS are designed to be as close as possible to 
the date of maximum luminosity. The few l ate UV spec tra we 
have in our sample come from IUE database (INES 2006). 



4. Training the model 

4. 1 . The training procedure 

The convergence process consists in minimizing a y 2 that per- 
mits the comparison of the full data set with the model of equa- 
tion Q] For each SN, the parameters are the normalization and 
coordinates along the principal components (xjO, a color and 
re-calibration parameters for spectra if any. The actual com- 
ponents (Mjfe) and the parameters of the color-law CL(A) also 
have to be estimated. This procedure requires a first guess for 
the model components (M*), for a first estimate of normaliza- 
tion, spectra re-calibration and color. We used the SALT model 
SED sequence for a SN with stretch = 1 for Mo and the dif- 
ference of SED sequence of a SN with stretch = 1.1 and the 
previous one for Mi (i.e. a linearized version of SALT model). 
Additional components where initiated with the orthogonal part 



of the SALT model SED sequence with respect to all previous 
components. 

We end up with more than 3000 parameters to fit, with ob- 
vious non-linearities, so that we used the Gauss-Newton proce- 
dure, which consists in: 

1 . Approximating locally the x 2 by a quadratic function of the 
parameters. 

2. Solving a large linear system to get an increment of the pa- 
rameters (5 Pi). 

3. Increment the parameters and iterate until the^ 2 decrement 
with respect to the previous iteration becomes negligible. 

First, the average model is estimated along with the color- 
law, calibration coefficients for spectra, and parameters of the 
SNe ( (x,), c ). When the system has converged, we add another 
component, and all the parameters are fitted again (components, 
color-law, SN parameters). The convergence algorithm is insen- 
sitive to the input set of components. 

4.2. Regularization 

There might be some degeneracy in part of the phase space for 
the given data set. For instance, if a phasexwavelength region 
is only covered by photometry and not spectroscopy, we do not 
have enough data to constrain the combinations of parameters 
that model spectral features, whereas we can still model a pho- 
tometric measurement, since the signal is integrated on a large 
spectral band. Adding a regularization term in the \ 2 solves this 
issue. If its contribution is low enough, it will not alter signifi- 
cantly the determination of parameters that are addressed by the 
data, while putting some limitation on the parameters that are 
not. We have chosen to minimize second derivatives with respect 
to phase and wavelength (once again, effective only when there 
is not enough data). The regularization term is the following : 



where M^ is the vector describing component k, D is the deriva- 
tive matrix and n a normalization that controls the weight of this 
regularization with respect to data. Since such a term introduces 
a bias in the estimator (departure from the maximum likelihood 
estimator), we have to quantify it in order to adjust the normal- 
ization n. For this purpose, we used a simulated dataset. This 
simulation helps us to define the resolution of the model. Each 
SN of the training sample was adjusted using the SALT model, 
then fake light-curves and spectra were computed by replacing 
each true measurement of the SN by the best fit value of the 
model. The training procedure applied to this data set gives a re- 
sult that is slightly biased due to the regularization term in the 
X 2 in the UV wavelength region. The weight of the regulariza- 
tion term (normalization «) was chosen so that the bias in K- 
corrections is smaller than 0.005 mag for all wavelength, which 
is significantly less than the statistical uncertainties. 

4.3. Model resolution 

The choice of the model resolution is imposed by the data set we 
have. We used 10 x 120 parameters for Mo (10 along the time 
axis and 120 for wavelength), in a phase range of [-20,-1-50] 
days and a spectral range of [2000, 9200] A. This gives a spectral 
resolution of order of 60 A which is sufficient for the modeling 
of SNe with broad lines due to the velocity of the ejecta. For Mi , 
we choose to use a lower resolution ( 10 x 60 parameters). The 
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time axis is remapped so that the time resolution at maximum 
and +20 days after maximum is a factor two better than at -20 
and +50 days (approximately 4.5 and 9 days respectively). As 
in SALT, we used only two free coefficients to model the color 
law CL(A) (third order polynomial, with two coefficients fixed so 
that CL(A B ) = and CL{A V ) = 0.4 log(10), see Eq.Q). Using this 
number of parameters, when the model is trained with the sim- 
ulated data set described above, we found that the limited reso- 
lution introduces a scatter in colors of only 0.01 mag standard 
deviation (it is a scatter rather than a systematic effect because 
of the varying epochs of photometric observations), which has 
a negligible on distances when compared to the intrinsic disper- 
sion of SNe la luminosities. 

5. Result of the training 

We decided to consider only two components for the current 
analysis since additional components are poorly constrained in 
most of the phase space and marginally significant in the re- 
gion of good data coverage. As the data sets improve so will the 
power to extract additional components. As a consequence, for 
each SN, we ended up with four parameters, a date of B-band 
maximum, a normalization, the parameter x\ and a color. The 
average value of x\ and its scale are arbitrary since we can mod- 
ify the components in consequence. We adopted < x\ >= and 
< x\ >= 1. 

Figure [2] shows the variation of the UBVRI light-curves as 
a function of parameter x\. We find that most of the variability 
can be described by a simple stretching of light-curves despite 
the fact that we did not force such behavior in the model. More 
quantitatively, the parameter x\ can be converted into a stretch 
factor^, whose actual value depends on the refe rence light-curve 
templa te used, here the one of SALT and of iGoldhaber et all 
d2001l) (ff-band ligh t curve template "Parab -18", G01); or into 
Am 15 dPhillipsll 19931) using the following transformations: 

s (SALT) = 0.98 + 0.091jq + 0.003x? - 0.00075x? 
s (G01) = 1.07 + 0.069jci - 0.015x? + 0.00067x? 
Amis = 1.09-0.161xi +0.013*? -0.00130^ 

Since there is not a perfect match of the non-linear stretch 
and A/7ii5 models with this one, those transformations (obtained 
with simulations) vary with the weight attributed to each phase 
(the scatter for stretch is about 0.02). 

We also notice the U - B color variation as a function of B- 
band light-curve broadening. The value of U - B for phase=0 
does not vary with X\ (see Fig. but when the flux is inte- 
grated in the phase range -10, +10 days, we find that U - B oc 
-0.2xs(SALT). This is about half the value obtained with SALT. 
However, we find a color law very close to the one obtained with 
SALT (see Fig. [3}, despite the fact that the supernova model is 
significantly different and the training set much larger. 

These results confirm the main findings of SALT. More inter- 
esting is the variability of spectra with the first component as dis- 
played in Figure |4] for three phases about maximum. It appears 
that the variability in U-band at maximum identified with light- 
curves can be attributed to a sharp variation of the spectrum for 
wavelengths lower than 3400 A. It is possible to identify such a 
feature thanks to the high redshift SNLS spectra. Indeed, the cal- 
ibration of ground-based spectroscopic observations in the near 

3 Of course, since the model is a linear combination of two compo- 
nents, we do not retrieve exactly the stretch model. 



ultraviolet may not be reliable, whereas UV spectra obtained by 
the International Ultraviolet Explorer have a very low signal to 
noise ratio for wavelengths larger than 3000 A. This feature can 
hardly be approximated by a broad-band color evolution. Hence 
we expect a net improvement of the accuracy of distance esti- 
mates in the UV range with respect to SALT or other equivalent 
methods which rely on a single spectral sequence. As an exam- 
ple of the modeling of the variability of spectral fea tures, fig ure[5] 
presen ts the variation of the R(Si II), as defined in lNugent et al.l 
(fl995l) . as a function of Am 15 retrieved fro m the model. It iscom- 
pared to a compilation of observations by iBenetti et al.l d2004l) . 
There is a good match in the Amis range of the training sample 
(Amis < 1.6). 

One must however evaluate carefully the statistical signifi- 
cance of the trends in the model. Statistical errors of this model 
rely on the weighting applied to spectra, and are sensitive to 
the errors assumed for photometric measurements which are not 
very secure for most nearby supernovae. Hence this accuracy of 
the modeling must be evaluated with distributions of residuals as 
described in the next paragraph. 

6. The remaining variability 

In order to assert the predictability of the model we need an 
independent data set that is not used in the training procedure. 
However, since we do not have a large number of measurements 
available, we resorted to a jackknife procedure: for each SN, 
we trained the model using all SNe from Table [2] but this one, 
and looked at residuals of the SN measurements to the retrieved 
model, fitting only the parameters concerning this particular SN, 
i.e. the date of maximum, (xo,i), and color. 

The residuals obtained by this method are a priori highly cor- 
related, and this correlation is difficult to estimate from first prin- 
ciples. However we would like to use this information to extract 
some intrinsic variability of the SNe, beyond the principal com- 
ponents of the model, in order to weight data according to this 
variability in addition to measurement errors. We resorted to the 
following simplification, using two kinds of residuals (or model 
errors) : i) Diagonal errors of the model are estimated from fits 
of light-curves with an independent normalization for each of 
them. The residuals obtained are of course still correlated, but 
we allow ourselves to treat them as independent, assuming that 
the correlation length (along time axis) is smaller than the data 
sampling for most light-curves Q. ii) /^-correction uncertainties, 
which can be estimated using the difference between the peak 
magnitude obtained from a single light-curve fit, and the one 
predicted by the model in the same filter, fitting all light-curves. 

6.1. Diagonal uncertainties 

We use the correlations between the estimates of the components 
given by the correlation matrix retrieved at the end of the train- 
ing procedure. However, we allow ourself to scale these statis- 
tical uncertainties on the model in p, A bins to account for the 
remaining variability. 

The model variance can then be defined as : 

V M odel(xi,p,A) = S(p,A)xV M ean{x\,p,A) 
Vmean{xx, P ,A) = H^VoHo + x?H[ViHi + 2x 1 HjC ,iHi 

where Ho(p, A) and Hi(p, A) are the vectors defined in Sec. I2.ll 
for components and 1, Vo, Vi and Co,i are the full variance and 

4 It is actually the purpose of the principal component analysis to 
extract all the correlations between observables for a given SN 
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covariance matrices of the components, and S(p,A) the scaling 
function. 

In each p, A bin, S (p, A) is evaluated so that 

1 [,/;■- xo^.Mq + ^H^.M!)] 2 _ 
nt-J 0-2 + xIS(p,A)Vmean{x\,P,A) 

where (fi) and (cr,-) are the measurements and their associated 
uncertainties. We evaluated separately these errors for light- 
curves and spectra, in order to take into account the correlated 
errors along the wavelength axis when dealing with photometric 
data. 

Photometric residuals of the jackknife-like procedure are 
shown figure [6] The data set are split in four rest-frame wave- 
length ranges, [3200,3900], [3900,5000], [5000,5700] and 
[5700, 7300] A, that roughly correspond to U, B, V, R and I 
bands respectively. The model uncertainties basically follow the 
statistical errors of data, with very good accuracy at peak bright- 
ness and a poor quality at early and late phases, especially in 
the [/-band. In the rise time region, the large errors are partly 
due to the limiting resolution of the model. Those errors are also 
displayed in figure|2]for |xi|=2. 

When fitting spectra, we imposed the values of the date of 
maximum, x\ and color obtained with the light-curves fit. The 
only remaining free parameters were those used to photomet- 
rically "re-calibrate" spectra. We found that model accuracy is 
poor at early phases and in the UV region (Fig. |4j». 

These estimates of the model errors can be accounted for 
when fitting the light-curves or spectra. It gives more reliable 
statistical errors for parameters (peak brightness, color, x\ and 
date of maximum) than when only statistical errors of measure- 
ments are considered. 

6.2. K-correction uncertainties 

A direct approach to access the quality of /^-corrections is to 
compare the observed peak magnitude of a light-curve in a given 
filter with the one predicted by the model using a fit of the other 
light-curves. Figure Q presents the differences between the ob- 
served and predicted magnitudes as a function of the effective 
rest-frame wavelength of the instrument response used. 

A more elaborated approach consists in modeling K- 
correction errors with a parametric function of wavelength 
which value vanishes for wavelength corresponding to the rest- 
frame B and V -bands (errors on B and V magnitudes at maxi- 
mum enter in the normalization and color evaluation). For each 
SN with enough light-curves, those /f-correction additional pa- 
rameters can be estimated and their standard deviation used to 
derive a model of /^-correction errors. Such a model is repre- 
sented by the solid line figure [7] and is given by the following 
formula: 

o- K (A)= 0.022 (j^f forA<A B 
= 0-018 (S) 2 A > A V 

Since the estimate of ctk is based on the fit of the normal- 
ization of light-curves, it measures a dispersion of colors aver- 
aged over the phase range defined by the data set. Clearly, the 
/if-correction errors are large in the UV range. Also, those er- 
rors must be added to the statistical errors on normalizations, 
but they do not account for the whole observed scatter. For in- 
stance, for An ^ 3600A, we find a dispersion of 0.04 magnitude 



(consistent with the one obtained in A06, Fig. 11.), but an uncer- 
tainty of only 0.022 magnitude has to be added to the statistical 
errors to match the observed dispersion. Diagonal uncertainties 
and /if-correction uncertainties are taken into account in the fits 
performed in the following section. 



7. Improving the distance estimates of distant SNe 

Despite the fact that our modeling is less accurate in the UV 
range, it is still very useful for distance estimates of high-redshift 
SNe (z > 0.8) for which, as in the case of SNLS, the rest-frame 
B and V-band observations often have a very poor signal to noise 
ratio. 

7. 1 . Light-curve fit of distant supernovae 

Figure [H shows the SN la SNLS-04D3gx at z = 0.91 fitted by 
the model. All four light-curves (g,r,i,z) are well described by 
the best fit model for which only four parameters were adjusted 
(date of maximum, normalization, color and x{). The^ 2 per de- 
gree of freedom (d.o.f.) of the fit is 0.76 (for 50 d.o.f.) when 
diagonal and /^-correction uncertainties of the model are con- 
sidered. Table Q] illustrates the gain in the accuracy of the color 
estimate for SNLS-04D3gx. 

bands /l ra j n color 

~Tz 3980 -0.220 ± 0.180(s) ± 0.033(d) ± 0.005(k) 

r,i,z 3250 -0. 147 + 0.05 1 (s) ± 0.026(d) ± 0.049(k) 

g,r,i,z 2520 -0.172 + 0.047(s) + 0.023(d) + 0.047 (k) 

Table 1 Error on the color estimate of SNLS-04D3gx as a func- 
tion of the number of light-curves included in the fit. The contri- 
butions to the error are measurement statistical errors (s), diago- 
nal model errors (d) and /^-correction errors (k). 

The total uncertainty on the color parameter is reduced by a 
factor 2.5 when g and r-band (rest-frame UV) light-curves are 
included in the fit. We see that the model uncertainties are large 
in this wavelength range and therefore must be propagated. 

7.2. Improving cosmological results 

The parameters retrieved from the light-curve fit can be used to 
estimate distances using the same procedure as described in A06. 
The distance estimator is a linear combination of m* B , x\ and c : 

I^b — m* B — M + a x X xi — B X c 

with m* B , x\ and c derived from the fit to the light curves, and 
ttxQ P an d the absolute magnitude M are parameters which are 
fitted by minimizing the residuals in the Hubble diagram. As in 
A06, we introduce an additional "intrinsic" dispersion (cr,„ f ) of 
SN absolute magnitudes to obtain a reduced x 1 of unity for the 
best fit set of parameters. 

We minimize the following functional form, which gives 
negligible biases to the estimates of a x and B: 

2 = y V T X S -M-5 log 1Q (</ £ (e,z)/10pc) 
x Zj v t C(X s ) V 



5 Note that the definition of a x differs from that of a in SALT 
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stands for the cosmological parameters that define the fitted 
model and di is the luminosity distance. C(X S ) is the covariance 
matrix of the parameters X s for which we have included in the 
variance of m* B the intrinsic dispersion and an error of the dis- 
tance modulus due to peculiar velocities, which we take to be 
300 km. s _1 . 

To estimate the systematic effect due to modeling of the 
SN la SED sequence, we fit the data set of A06, which con- 
sists in 44 nearby SNe la and 71 SNLS SNe. Whereas we used 
improved photometry for the training of the model, we consider 
here exactly the same data set as in A06, in order to ease the 
comparison with the results obtained with SALT. 

Thanks to the evaluation of the model errors, we can safely 
use all available light-curves in the fit. Especially the r-band 
data at redshifts greater than 0.8 (effective rest-frame wavelength 
lower than 3440 A) are very useful to constrain the color of the 
supernova. With this additional information, the uncertainty in 
distance moduli is significantly reduced at high redshift, yielding 
a better resolution on cosmological parameters. For flat ACDM 
cosmology, we obtain : 

Q M = 0.240 + 0.033 
a x = 0.13 + 0.013 
P = 1.77+0.16 

with cr im = 0.12 (which is smaller than the value of 0.13 in 
A06 because we have considered uncertainties in the model). 
The RMS of the residuals around the best fit Hubble relation is 
0.161 mag (compared to 0.20 in A06, see Fig. [9]). 

The uncertainty on fl^ is improved by 10% with resp ect to 
the A06 (Q M = 0.263 + 0.037, table 3 of lAstier et ai1l2006l) anal- 
ysis. We find a difference of 0.023 on Q.m which is consistent 
with the assumed systematic error due to modeling of 0.02 in 
that paper. However the supernova SNLS-03Dlcm at a redshift 
of 0.87 now appears as a significant outlier as the distance mod- 
ulus resolution improved^. This object has a spectrum with low 
signal to noise ratio and was classified as probable la. Discarding 
this object from the analysis gives fl^ = 0.246 + 0.032 and a 
standard deviation of residuals of 0. 154 magnitude. 

Since the current model significantly improves distance es- 
timates at high redshifts, we obtain a greater improvement on 
the estimation of the equation of state of dark energy. As an ex- 
ample, we may us e the figure of merit proposed by the Dark 
Energy Task Force dAlbrecht et al.l2006b . which is inversely pro- 
portional to the area of the 95% confidence level contour in the 
plane (w„, w p ), where the following parametrization is consid- 
ered for the equation of state of dark energy : 

w = w p + (a p - a)w a 

a being the scale factor, and a p a reference scale factor chosen 
so that the estimates of w p and w a are uncorrected for a given 
expe riment. Using th e baryo n acoustic oscillations constraints 
from [Eisenstein et all d2005l) (Eq. 4), a p = 0.851, and we im- 
prove this figure of merit by 35% with respect to the analysis 
using SALT (for which a p is slightly different). 

6 SNLS-03Dlcm is 0.6 ± 0.2 mag dimmer than expected for the best 
fit cosmology, which corresponds to a 3 <x deviation when including the 
intrinsic dispersion. This has a 27% probability to occur at least once by 
chance for our sample of 1 15 SNe, if SNe distances are scattered about 
the Hubble law following a Gaussian distribution. 



Light-curves can be fitted with all the models from the jack- 
knife procedure, so that we can derive as many estimates of Q. M 
as there are SNe in the training sample. We found that the RMS 
of this distribution is 0.003, so that we expect a deviation from 
a model trained with an infinite number of SNe of the order of 
0.0015 on Q.m (for a flat ACDM cosmology), a value that is neg- 
ligible compared to the other sources of systematic errors in such 
an analysis. 

8. Other applications 

The proposed model provides a tool for spectroscopic and pho- 
tometric identification, and photometric redshift determination. 
A detailed analysis of the purity and efficiency of a photometric 
identification tool with respect to SNe lb, Ic and II is beyond the 
scope of this paper. 

8.1. Spectroscopic identification 

The model proposed allows simultaneous fits of light-curves and 
spectra (with additional galaxy templates to evaluate the host 
contamination as mentioned in Sec. I3.2I ). This can turn out to 
be very useful for identification of Type la supernovae. Indeed, 
some Type Ic SNe can present light-curves and spectra that look 
qualitatively like SNe la, and in the most extreme cases, pho- 
tometric and spectroscopic identifications taken separately may 
fail to tag this object correctly. A detailed analysis is beyond the 
scope of this paper. We will show only two examples. 

First, figure \TU\ shows the observed light-curves and spec- 
trum of SNLS-03D4ag at a redshift of 0.285 along with the re- 
sult of the simultaneous fit. Four "re-calibration" parameters for 
the spectrum were considered (since there are four light-curves). 
The x 1 per degree of freedom for the light-curves and the spec- 
trum are respectively 0.53 (for 28 d.o.f.) and 0.63 (for 1770 
d.o.f.), taking into account the model errors Q), so that this SN 
can be safely considered as a typical normal SN la. 

Now if we consider an SN Ic, for instance SNLS-03D4aa at 
a redshift of 0. 166, the model gives a very bad fit of the data (re- 
duced x 2 of 4.6 (for 10 d.o.f.) and 1.6 (for 1807 d.o.f) for light- 
curves and spectra respectively), allowing for a clear rejection of 
this event. 

8.2. Photometric redshifts 

For photometric redshift determination, one can compare the 
redshift estimate based on a simultaneous fit of all light-curves 
of a given supernova and the much more precise spectroscopic 
redshift derived from the spectroscopic observation of the host of 
the object. Using all available light-curves (up to four for SNLS) 
allows us to relax assumptions on the light curve shape parame- 
ter (jc i ) and the color, without any use of the absolute luminosity, 
so that we do not require any prior on the cosmological parame- 
ters. 

We have applied this method to all SNe la from Tablef2]with 
at least 3 light-curves, using for each SN the model obtained 
without this SN (always to avoid over- training). A Gaussian 
prior for x\ was assumed based on the statistics from the training 
procedure, x\ — ± 1. No prior was applied to color. 

The resulting photometric redshift is compared to the spec- 
troscopic one on Figure QT| The RMS of the distribution of 
Az/(1 + z) is 0.02 with no significant bias. A Gaussian fit of 

7 With the model errors, the average x 1 P er degree of freedom for all 
the SNe of the training sample is one by definition 
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the distribution of Az/(1 +z) gives cr = 0.01. This is an accuracy 
of the same order of magnitude as the one that can be derived 
from a fit of the SN spectrum alone. The only low-redshift SN 
for which the photometric redshift is off by more than 0.1 is 
SN1999cl which is a highly extincted supernova. 

9. Conclusion 

We have proposed a new empirical model of Type la super- 
novae spectro-photometric evolution with time, based on ob- 
served nearby and distant supernovae light-curves and spectra. 
The method uses available spectra, regardless of their wave- 
length range or calibration accuracy, since they are re-calibrated 
using the photometric information in the training process. This 
model provides an average spectral sequence of Type la super- 
novae and their principal variability components^. 

Thanks to an evaluation of the modeling errors (for photo- 
metric points, spectra, and broad-band colors), on can safely use 
most of the information on a given supernova for comparison 
with the model. This is very helpful for distance estimates, but 
also photometric redshift evaluation and SN identification. 

Applied to the supernovae sample of A06, we improved the 
distance estimates, especially at redshifts larger than 0.8, thanks 
to the modeling of the ultraviolet emission and the propagation 
of model errors. When a flat ACDM cosmology is fitted to the 
Hubble diagram, the gain in statistical resolution on fl^ is of 
10% compared to the previous analysis using SALT; it is 35% 
better on the two parameter constraint of the equation of state of 
dark energy. 

The model accuracy is currently limited by the size and qual- 
ity of the current data sample. A larger sample would be needed 
to look for another variability component. 

This empirical modeling technique is well suited for 
the analysis of very large samples of supernovae photo- 
metric data such as the ones ex pected fr om future proje cts 
like DUNE dRefregier et al.ll2006h . JDEM dAlderinglliool or 
LSST 0. Indeed, with large statistics of high precision photo- 
metric data spread on a large redshift range, we may expect to 
deconvolve the SED sequence of the average supernova and its 
principal variations. Such an analysis would permit one to look 
for variability correlated with redshift, and hence provide a non- 
parametric approach to test for supernova evolution. 
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Fig. 3 The color law c x CL(A) as a function of wavelength for 
a value of c of 0.1 (solid line). The dashed curve represents the 
extinc tion with respect to B band, {A a —Ab), from lCardelli et all 
(119891) with R v = 3.1 and E(B - V) = 0.1, and the dotted line 
is the color law obtained with SALT (very close to the result 
obtained here). 
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Fig. 2 The U*UBVRI template light curves obtained after the 
training phase for values of x\ of -2, 0, 2 (corresponding to 
stretches of 0.8, 1.0, and 1.2; dark to light curves) and null By 
color excess. U* is a synthetic top hat filter in the range 2500- 
3500 A. The shaded areas correspond to the one standard devia- 
tion estimate as described in section I6TT1 
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Fig. 5 T he solid curve presents the relation of 
R(Si II) (iNugent et all 1 19951) as a functio n of Am 15 along 
with measurements from lBenetti et alJ (12004 . 
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Fig. 4 Spectra at -10, 0, +15 and +35 days about B-band max- 
imum for values of x\ of -2, 0, 2 (corresponding to stretches 
of 0.8, 1.0, 1.2; light gray, black curve, dark gray) and null By 
color excess. The shaded areas correspond to the one standard 
deviation estimate as described in section loTTI The dashed curves 
represents the spectra of N02 (version 1.2). All spectra are arbi- 
trarily normalized. 
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Fig. 6 Estimated standard deviation of model photometric errors 
as a function of phase, for several rest-frame wavelength ranges 
roughly corresponding from top to bottom to U, B, V, R and 
/-bands. Those model errors were evaluated from the scatter of 
residuals to the single light-curve fit. 
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Fig. 7 Difference between observed peak magnitude in each band 
of each SN from tableland the model prediction as a function 
of the rest-frame effective wavelength of the filter used (gray tri- 
angles : SNLS SNe, gray squares : nearby SNe). The large black 
symbols represent the estimated dispersion in each wavelength 
bin (triangles for SNLS, and squares for nearby SNe). The large 
circles show the average difference in each wavelength bin for all 
SNe and the solid curve is a polynomial fit to the dispersion used 
as an estimate of the /^-correction scatter. Since uncertainties on 
B and V magnitudes at maximum enter in the normalization and 
color evaluation of the model, /f-correction uncertainties are set 
to zero for B and V- band wavelengths. 
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Fig. 8 Observed light-curves points of the SN la SNLS-04D3gx 
at z=0.91 along with the light-curves derived from the model 
(solid line, trained without this SN). The dashed lines represent 
the 1 cr uncertainties of the model (both uncorrelated and K- 
correction errors). 
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Fig. 9 Residuals to the distance redshift relation for the best fit 
ACDM cosmology for SNLS and nearby SNe la. 
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spec-z 

Fig. 1 1 Photometric vs spectroscopic redshifts for the training 
sample of SNe. The dashed lines are z + 0.1(1 + z) functions. 
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Fig. 10 Top: Observed light-curves points of the SN la SNLS- 
03D4ag at z=0.285 along with the light-curves derived from the 
model (solid line). Bottom: Observed VLT spectrum at a phase 
of -8 days along with the model spectrum. The dashed curve 
is the expected model from the light-curve fit (up to a global 
normalization) and the solid curve the one including four re- 
calibration parameters. 
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Table 2. The training sample of Type la supernovae light curves and spectra 



Name 


z d 


Bands' 5 


Number of spectra; phases c 


1981B 


0.006 


UBV (B83) 


6 ; < p < 35 (Br83) 


1986G 


0.002 


BV (P87) 


6; -3 < p < 18 (AC) 


1989B 


0.002 


UBVR (W94) 


18 ; -6<p< 19 (AC) 


1990af 


0.051 


BV (H96) 





1990N 


0.003 


UBVR (L91,L98) 


15; -14 < p < 38 (AC,IUE,CFA) 


1991T 


0.006 


UBVR (F92,L98,A04,K04b) 


36 ; -13 < p < 258 (AC,CFA) 


1992 A 


0.006 


UBVR (S92.A04) 


19 ; -5 < p < 27 (CFA.IUE) 


1992al 


0.015 


BVR (H96) 





1992bc 


0.020 


BVR (H96) 





1992bo 


0.019 


BVR (H96) 





19930 


0.052 


BV (H96) 





1994ae 


0.004 


BVR (R99,A04) 





1994D 


0.001 


UBVR (R95,P96,M96,A04) 


21 ; -11 < p<76(AC) 


1994S 


0.015 


BVR (R99) 


1 ; p = 22 (AC) 


1995ac 


0.050 


BVR (R99) 





1995al 


0.005 


BVR (R99) 





1995bd 


0.016 


BVR (R99) 





1995D 


0.007 


BVR (R99,A04) 





1995E 


0.012 


BVR (R99) 





1996W 


0.036 


BVR (R99) 





1996bo 


0.017 


UBVR (R99,A04) 





1996X 


0.007 


BVR (R99) 


12 ; -1 <p< 57 (AC) 


1997bp 


0.008 


UBVR (J05) 





1997bq 


0.009 


UBVR (J05) 





1997do 


0.010 


UBVR (J05) 





1997E 


0.013 


UBVR (J05) 





1998ab 


0.027 


UBVR (J05) 





1998aq 


0.004 


UBVR (R05) 


26 ; -9 < p < 90 (CFA) 


1998bu 


0.003 


UBVR (S99) 


66 ; -3 < p < 56 (CFA) 


1998dh 


0.009 


UBVR (J05) 





1998ef 


0.018 


UBVR (J05) 





1998es 


0.011 


UBVR (J05) 





1999aa 


0.014 


UBVR (J05) 


9 ; -11 <p <40 (G04) 


1999ac 


0.009 


UBVR (J05) 





1999aw 


0.038 


BVR (St02) 





1999cl 


0.008 


UBVR (J05,K00) 





1999dk 


0.015 


UBVR (K01, ACM) 





1999dq 


0.014 


UBVR (J05) 





1999ee 


0.011 


UBVR (S02) 


11; < p < 43 (H02) 


1999ek 


0.018 


UBVR (J05,K04b) 





1999gp 


0.027 


UBVR (J05,K01) 





2000E 


0.005 


UBVR (V03) 





2000fa 


0.021 


UBVR (J05) 





200 lba 


0.029 


BV (K04a) 





200 lbt 


0.014 


BVR (K04b) 





200 lcz 


0.016 


UBVR (K04b) 





200 lei 


0.004 


UBVR (K03) 


2; p = 24,33 (W03) 


2001V 


0.015 


BVR (Vi03) 





2002bo 


0.004 


UBVR (Z03,B04,K04b) 


12 ; -14 < p < 38 (B04) 


2002er 


0.009 


UBVR (P04) 





2003du 


0.006 


UBVR (A05) 


4 ; -11 < p < 37 (A05) 


2004fu 


0.009 


UBVR (T06) 





SNLS-03Dlau 


0.504 


griz 





SNLS-03Dlaw 


0.582 


griz 
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Table 2 (cont'd) 



Name z a Bands' 5 Number of spectra; phases 1 



c\tt c n arn qv 

OiN Lo-UjIJ 1 aX 


u.^-yo 


griz 


i . 
i , 


P = 


—j 




a o^c 
0.803 


griz 


1 : 


P = 


—o 


olNLo-UJJJ ICO 


u.o/y 


griz 


i . 


P = 


i 

—4 


CMT C Pil T^i 1 «... 


a o^o 
0.805 


griz 


1 


; p z 


t 

= 1 




a 1 1 1 
0.331 


griz 


i . 
1 ; 


p = 


—5 


qmt c mn 1 fi 


U.DOO 


griz 


i 
i 


; p z 


- 1 




0.800 


griz 









oINLo-U.3L)4ag 


a ooc 


griz 


i 

1 ; 


p = 


o 

-o 


CMT C HI T~"\ A 


0.633 


griz 


1 


; p z 


- 6 




a OTA 

O.Z/U 


griz 


i . 
1 : 


p = 


—y 


CMT c mr\/i n . r 
5JNL5-Uj1J4cx 


a n/in 


griz 








CMT C flQn/l,... 

oiNLo-U-jD4cy 


o.y// 


griz 




A 




cmt c mn/i/.r. 


O.oya 


griz 




n 





cmt c mn/i/iu 


O.Oz / 


griz 








CMT C (\1T\AA,, 

oJNLo-UjJJ4ay 


n AA/i 
0.6U4 


griz 


1 


; P = 


A 

= 4 


CMT C P\1T\A£A 


a 7n i 

o. /yi 


griz 


i 

1 ; 


P = 


t 

— 1 


CMT C AQr\/1/>l 

oJNLo-UJlMgl 


0.5/1 


griz 








CMT C f\AT\ \ n~ 

oJNLc>-U4Dlag 


0.55 / 


griz 


1 


; p z 


= 4 


CMT C C\ATW An 


O.zl 1 


griz 


i 

1 ; 


p = 


t 

— 1 


CMT C n /1 1 £F 


a otn 
O.80U 


griz 


i 
1 


; p z 


A 

= 4 


CMT C r\AT\1UA 


U.3uy 


griz 


i 
t 


; p z 


- U 


CMT C r\/1T~\1l-.,r 


n caa 
0.56U 


griz 


1 


; p z 




cmt c n/ir>iii., 
jJN Lo-U41J 1 ny 


a qca 


griz 


i . 


p = 


— 1 


cmt c n/n~\i:. r 


o.yys 


griz 


i 
1 


; p z 


= 1 


cmt c r\/iT~\ 1 
5JNL5-U4D1KJ 


a CO/1 
0.384 


griz 


i . 
1 : 


p = 


A 

—4 


CMT C Pi A 1 


o. /ys 


griz 


i . 
1 ; 


p = 


t 

— 1 


CMT C fT/1T"\1^-.U 


n caa 


griz 




u 




CMT C A/1T~\1^,.,7 


o.y id 


griz 


i 
1 


; p -- 


c 

- J 


C\TT C f\A p» 1 

vMNJ^o-U41J IpC 


U. / jo 


griz 




u 




CMT C Pi A T^\l n A 

3lNLiS-U4Dlpa 


n aco 

o.ysz 


griz 








CMT C A/1 1 

MNLc>-U4lJlpg 


U.J i3 


griz 


i . 


p = 


1 

— 1 


CMT C A,lfM ^„ 

oINLo-U4Dlpp 


/me 
0. /35 


griz 




A 

u 




CMT C A/IHUo 


a con 

u. joy 


griz 


t . 


p = 


—L 


cmt c ri/ir\i 


n TA/i 
0. /U4 


griz 




A 




CMT C H/1 T^ 1 nl* 


U.DOJ 


griz 




A 




CMT C Pi A T^O^w^ 

o IN Lu b - U4D Z ac 


a "2/1 o 


griz 




A 




CMT C C\A r^o^f 


u.ooy 


griz 


i 
t 


; p z 


o 

- <S 


cmt c r\/i r»ofn 
5JNLo-U4DZrp 


n yl 1 c 

0.413 


griz 


1 


; p z 


t 

= 1 


CMT C fT/1 T^Of^ 

olNLo-U4lJZrs 


U.3J / 


griz 


i 
t 


; p = 


i 

- 1 


CMT C A/1T~\O^U 

oINLo-U4DZgb 


n /i in 
0.43U 


griz 




A 




.MNL>3-U4DZgC 


n co i 


griz 


i . 


P = 


i 

—4 


CMT C r\y1T~\0^r. 

o JN L 5 - U4D 2gp 


0. /U/ 


griz 


1 


; p = 


= z 


CMT C Pi a r^oi^.- 


n T/i /i 
0. /44 


griz 




A 




CMT C A/irM/,^ 


a Aon 
O.DzU 


griz 




A 




SNLS-04D3dd 


1.010 


griz 


1 


; p = 


- 2 


SNLS-04D3df 


0.470 


griz 









SNLS-04D3do 


0.610 


griz 









SNLS-04D3ez 


0.263 


griz 









SNLS-04D3fk 


0.358 


griz 









SNLS-04D3fq 


0.730 


griz 


1 


; p z 


= 1 


SNLS-04D3hn 


0.552 


griz 









SNLS-04D3kr 


0.337 


griz 


1 


; p = 


= 4 


SNLS-04D3mk 


0.813 


griz 









SNLS-04D3ml 


0.950 


griz 
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Table 2 (cont'd) 





7 a 


Rands b 


Number of spectra; phases 


SNLS-04D3nh 


0.340 


griz 


t 
1 


; p = 3 


SNLS-04D3ny 


0.810 


griz 


1 


; p= 1 


SNLS-04D3oe 


0.756 


griz 




u 


SNLS-04D4an 


0.613 


griz 




u 


SNLS-04D4bq 


0.550 


griz 


L , 


p = z, 4 


SNLS-04D4dm 


0.811 


griz 




U 


SNLS-04D4fx 


0.629 


griz 


1 ; 


p = -8 


SNLS-04D4gg 


0.424 


griz 




u 


SNLS-04D4hu 


0.703 


griz 







SNLS-04D4ic 


0.680 


griz 


i 
t 


> p - 2 


SNLS-04D4ii 


0.866 


griz 




u 


SNLS-04D4im 


0.751 


griz 




u 


SNLS-04D4in 


0.516 


griz 







SNLS-04D4jr 


0.482 


griz 


1 ; 


p = -6 


SNLS-05Dlck 


0.617 


griz 







SNLS-05D2ab 


0.320 


griz 




u 


SNLS-05D2ac 


0.494 


griz 







SNLS-05D2ah 


0.180 


griz 




u 


SNLS-05D2ay 


0.921 


griz 







SNLS-05D2bt 


0.679 


griz 







SNLS-05D2bv 


0.474 


griz 




u 


SNLS-05D2bw 


0.921 


griz 




{) 


SNLS-05D2ci 


0.631 


griz 







SNLS-05D2cl 


0.839 


griz 




u 


SNLS-05D2ct 


0.734 


griz 







SNLS-05D2dw 


0.417 


griz 







SNLS-05D2dy 


0.498 


griz 







SNLS-05D2eb 


0.639 


griz 




r\ 



SNLS-05D2ec 


0.526 


griz 




r\ 



SNLS-05D2fq 


0.735 


griz 




r\ 
K) 


SNLS-05D2hc 


0.360 


griz 







SNLS-05D2he 


0.608 


griz 




u 


SNLS-05D3ax 


0.643 


griz 




u 


SNLS-05D3cf 


0.419 


griz 




r\ 
U 


SNLS-05D3ci 


0.515 


griz 




u 


SNLS-05D3cx 


0.805 


griz 




u 


SNLS-05D3dd 


0.480 


griz 




u 


SNLS-05D3gp 


0.580 


griz 




u 


SNLS-05D3gv 


0.715 


griz 




u 


SNLS-05D3ha 


0.808 


griz 




u 


SNLS-05D3hq 


0.338 


griz 




u 


SNLS-05D3hs 


0.664 


griz 




r\ 



SNLS-05D3ht 


0.900 


griz 




r\ 
K) 


SNLS-05D3jb 


0.740 


griz 




u 


SNLS-05D3jh 


0.718 


griz 







SNLS-05D3jk 


0.736 


griz 







SNLS-05D3jq 


0.579 


griz 







SNLS-05D3jr 


0.370 


griz 







SNLS-05D3kp 


0.850 


griz 







SNLS-05D3kt 


0.648 


griz 







SNLS-05D3kx 


0.219 


griz 







SNLS-05D31b 


0.647 


griz 


1 ; 


p = -l 


SNLS-05D3mq 


0.240 


griz 







SNLS-05D4af 


0.502 


griz 
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Table 2 (cont'd) 







; a 


Rands b 


Numhpr of snpptra" nhasps c 


SNLS-05D4ag 





636 


griz 





SNLS-05D4av 





543 


griz 





SNLS-05D4be 





540 


griz 





SNLS-05D4bj 





704 


griz 





SNLS-05D4ca 





608 


griz 





SNLS-05D4cq 





783 


griz 





SNLS-05D4cs 





735 


griz 





SNLS-05D4cw 





375 


griz 





SNLS-05D4dt 





407 


griz 





SNLS-05D4dw 





853 


griz 





SNLS-05D4dy 





810 


griz 






"Heliocentric redshift. 



b Phot ometry 
P87: 



Referen ces : B83: iButa & Turned 
ap <fl987l) . W9 4: IWells et al.1 
H96: IHamuv et alT dl996al) . L91: iLeibundgut et al] 
L98: iLiraetalJ dl998h. F92: iFilippenko et al.l 
A04: lAltavilla et alj d2004l). K04b: iKrisciunas et all 




198 
199 
dl991h 
d!992h 
d2004bh 



S92:ISuntzeffldl992h. R99:lRiess et al.ldl999l) R95:lRichmond etal 



GUI 



Meikleetal.l~{ l996) 



, P96: iPatat et alj dl996|). M 96: 
J05:lJha et al"] J2006h. R05: iRiess et al] d2005T) , S9 9: ISuntzefF et al 
1999^ St02: IStrolger et al] d2002T). K00: IKrisciunas e~ 



iKrisc 

20001), K01: I Krisciunas et alj d200lh , S02: Stritzinger et al 
20021) V03: Ivalentini et al] d2003h. K 04a: IKrisciunas e~ 
2004id) . K03 : IKrisciunas et al J d2003h . Vi03: IVinko et al 
Z03: IZapata et al] d2003h . B 04: iBenetti et al] d2004j) 
d2004h . A05: lAnupama et ai] 120051) 



2003) 

P04: iPignata et al 



T06: lTsvetkovf 72006) 



c Only spectra with negligible hos t-contaminati on are 

listed . Spectroscop i c Ref erences : Br8 3 : iBranchetal] dl983f) 
AC: " 



Barbon et al 



2006), G04 



d!999h. IUE : 



Garavini et al] 



d2006|) , C FA: iBlondir 
H02: IHamuv et al. 



2002b, W03: IWangetal l d2003al) . B04: IBenetti et al] d2004l) 



A05: lAnupama et al] d2005n 



